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Abstract 

We prove that the eigenvalues of a certain highly non-self-adjoint 
operator that arises in fluid mechanics correspond, up to scaling by 
a positive constant, to those of a self-adjoint operator with compact 
resolvent; hence there are inflnitely many real eigenvalues which accu- 
mulate only at ±00. We use this result to determine the asymptotic 
distribution of the eigenvalues and to compute some of the eigenvalues 
numerically. We compare these to earlier calculations in [T], [2] and 
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1 Introduction 

In a recent paper [3j, we showed that the spectrum of the highly non-self- 
adjoint operator —iH is real, where H is the closure of the operator Hq on 
L^(— 7r,7r) defined by 

W)W=4(.inm|)+| (1) 

for any fixed s G (0,2) and all / G Dom (ifo) = Cpgr([— tt, tt]). Boulton, 
Levitin and Marietta subsequently proved in a recent paper that a wider 
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class of operators possess only real eigenvalues. However, they did not prove 
that any eigenvalues exist for these operators, nor that their spectra are real. 
The results obtained in this paper for the original operator ([T]) are much 
more detailed than those presented in [H |5]. 

The operator H was first studied by Benilov, O'Brien and Sazonov, who 
argued in [T] that the equation 

% = Hf (2) 

approximates the evolution of a liquid film inside a rotating horizontal cylin- 
der. They also made several conjectures, based on non-rigorous numerical 
analysis, including that the spectrum of H is purely imaginary and consists 
of eigenvalues which accumulate at ±zcx). 

Davies showed in [2] that —iH has compact resolvent by considering the 
unitarily equivalent operator A on /^(Z) defined by 

{Av)n = ^n{n - l)vn-i - ^n{n + l)vn+i + nVn (3) 

for all V G Dom(A) = [v E /^(Z) : Av E /^(Z)}. Here A = J^-\-iH)T, 
where JF : 7r,7r) — > /^(Z) is the Fourier transform. If JF/ = v then 

{vn)nez are the Fourier coefficients of /. This result was achieved by obtaining 
sharp bounds on the rate of decay of eigenvectors and resolvent kernels, and 
by determining the precise domains of the operators involved. He also showed 
that 

A = A_ © © A+, (4) 

where A_ and A^ are the restrictions of A to /^(Z_) and /^(Z+) respectively, 
and that A_ is unitarily equivalent to — Since the resolvent is compact 
and the adjoint has the same eigenvalues, the spectrum of —iH consists 
entirely of eigenvalues. 

As previously mentioned, we proved in [4J that these eigenvalues, if they 
exist, must all be real. Eigenvalues of H or —iH have been calculated nu- 
merically in [H [21 E] 5 but until now it has not been proven rigorously that 
any non-zero eigenvalues exist. 

In this paper we prove rigorously that —iH has infinitely many eigenvalues 
which accumulate at ±oo (Corollary 14.51) . Our approach is to show that the 
eigenvalues of A^ correspond, up to scaling by a positive constant, to those 
of a self-adjoint operator with compact resolvent (Corollary \2.2\ Theorem 
13.21 Theorem 14.21) . By analysing the self-adjoint operator, we determine the 
asymptotic distribution of the eigenvalues (Theorem 1 5. 6 1) . It was argued in [3] 
that the distribution of the eigenvalues, if they exist, should be quadratic, but 
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no rigorous bounds were given. We prove rigorously that A„ ~ evr^n^/^"^ for 
some constant j3 which we determine. We also perform numerical calculations 
of eigenvalues, which we compare to those given in [1], [2] and [3] (Section 
15.41) . Moreover, our calculated values are rigorous upper bounds on the 
true values of the eigenvalues, insofar as the computed eigenvalues of regular 
Sturm-Liouville problems, which are known to be computationally stable, 
can be said to be rigorous. This gives us some idea of the accuracy of the 
previous calculations. 

The correspondence of the eigenvalues of iH to those of a self-adjoint 
operator Q might lead us to believe that iH is similar to Q in the sense that 
there exists a bounded linear operator S with bounded inverse such that 
SDom (Q) = Dom(if) and iH = SQS~^. However, it has recently been 
proven that iH is not similar in this sense to any self-adjoint operator P, 
Proof of Theorem 5.1]. 

2 Correspondence of eigenvalues to those of 
a Sturm-Liouville problem 

We have already shown in [1] that if A is an eigenvalue of the operator 
defined on its natural maximal domain by 

then II = 2X/e is an eigenvalue of the Sturm-Liouville problem 

— {pu'Y = fiwu, (5) 

where 

p{x) = (l-a:)^+l/^(x + l)l-l/^ (6) 
w{x) = x-\l - xY^'ix + ly'/' (7) 

and u G C°°([0, 1]) with u{0) = 0. Moreover, the solution of ([5]) satisfying 
these conditions is 

oo 

«(^) = E ^nX"", (8) 
n=l 

where f„ is the solution of the recurrence relation 

Tl — A 

n{n - l)vn-i - n{n + l)vn+i + 2— — Vn = (9) 

satisfying the initial conditions vi = 1, f2 = (1 — A)/£. 
We now show the converse: 
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Theorem 2.1 If n is an eigenvalue of the Sturm- Liouville problem then 
X = efi/2 is an eigenvalue of A+. 

Proof If (f„) is the solution of the recurrence relation ([9]) satisfying the 
stated initial conditions and u is defined by on (0, 1), then m is a non-zero 
solution of 1^. Equation is equivalent to 

„ fl + l/e l-l/e\ , a 
u" + '— + '— u' - ^. -u = 0, 10 

so we see that a second linearly independent solution is U\ = au{z) log 2; + 
J2'^=o^nZ"', with 60 7^ 0. Suppose that /i is an eigenvalue of the Sturm- 
Liouville problem and y is a corresponding eigenvector. We proved in [1] that 
/i G R. Now y is a non-zero solution of (JSj) in (0, 1) such that \imx-,o+ y{x) = 
and lima;^i_ y{x) is finite. Since the space of solutions of ([5]) is two- 
dimensional, y = au + Pui for some a, P E C. Considering the end-point 
a; = 0, we see that we must have y = au. Without loss of generality, we may 
assume y = u. Hence u{x) converges to a finite limit as x — 1— . Suppose 
that A is not an eigenvalue of A^. Davies showed in [2] that, for A G R, (E]) 
has two linearly independent solutions 0, ip such that (j)n > n^^'^~^ > for 
all sufficiently large n and \il)n\ ~ n~^/^~^ as n ^ 00. The space of solutions 
of © is two-dimensional so f„ = a0„ -|- bipn, and a 7^ since ip ^ and 
V ^ /^(Z+). Without loss of generality a > 1. Hence there exists > such 
that Vn > for all n > N. For x G (0, 1), 

N-l 00 

uix) > ^ K - r2-i)a;" + ^ n-^x" 

n=l n=l 
iV-1 

= 5Z K - - log(l - x) 

n=l 
00 

as X 1 — . This is a contradiction, so A is an eigenvalue of A^. | 

Corollary 2.2 A is an eigenvalue of A^ if and only if fi = 2X/e is an 
eigenvalue of the Sturm- Liouville problem 



3 Self-adjointness 

We now show that the operator corresponding to the Sturm-Liouville problem 
is essentially self-adjoint on a suitable domain. Equation (E]) can be written 
as 

Lu = fiu (11) 
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where L is an operator on L^((0, l),w{x)(\x) defined by 

Lf^-w-\pr)' (12) 

on Dom(L) = Co-([0,l]) = {/ G C~([0,1]) : /(O) = 0} C L2((0, 1), w(a;)dx). 
We define ^ 

{f,9)w^ f{x)g{x)w{x)dx 

•J 

and ll/IL = for all /, g E L'\iO,l),wix)dx). 

We also consider Lc, the restriction of L to C^(0, 1), which is the space 
of smooth, compactly supported functions on (0, 1). 

Proposition 3.1 The adjoints L* of Lc and L* of L are closed extensions 
of L, which is symmetric, L* C L*, and the following are equivalent: 

(a) II E C is an eigenvalue of the Sturm-Liouville operator L and u e 
Dom (L) is a corresponding eigenvector; 

(h) 11 E C is an eigenvalue of the Sturm-Liouville operator L and u e 
Dom {l^ is a corresponding eigenvector; 

(c) fi G C is an eigenvalue of the operator L* and u e Dom (L*) is a 
corresponding eigenvector. 

Moreover, if < e < 1 then statements (a) - (d) are equivalent to 

(d) II & C is an eigenvalue of the operator L* and u G Dom (L*) is a 
corresponding eigenvector. 

Proof We first show that L is symmetric. For all f,gE Dom (L) we have 
{Lf,g)^ = - !\pf')\x)W)^x 

J 

= -p(l)/'(l)^ + p(0)/'(0)^+ C f\x)p{x)^dx 

Jo 



f{l)p{l)g'{l)-f{0)p{0)g'{0)- I f{x){pg'y{x)dx 
if, Lg)^ 



since p(l) = /(O) = ^(0) = 0. 

Now it is clear that Lc is also symmetric and 



Lc^LCL* C L*, 



the last two being closed. 

We now prove the equivalence of statements (a)-(c): 
(a) =^ (b) =^ (c): Immediate. 
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(c) (a): For all (f) G C^{0, 1) we have 



u{x){p(f)'y{x)dx 



■u, Lr 



{L^)*u,(t) 



—fi / u{x)(f){x)w{x)dx 



so {pu')' = —fiwu when we consider u as an element of the space T>'{0, 1) of 
distributions on the test-function space C^(0, 1). Since u G L^((0, 1), w{x)dx), 
u G L2(0, 1) C Li(0, 1). Also ^1/2 g (3(^(^0^ 1])^ gQ g 1) for any 
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5 G (0,1). Therefore pw' G 1) for any such 5 and hence pw' has a 

representation which is continuous on (0, 1] given by 

rl 



{pu')'{y)dy + c = jj, w{y)u{y)dy + c 



{pu') (x) 



for some constant c and all x G (0, 1]. Since p is continuous on [0, 1] and p > 
on [0,1), u' is continuous on (0,1), i.e. u is continuously differentiable on 
(0, 1). It now follows from the above equation that pu' is in fact continuously 
differentiable on (0, 1). Since p is continuously differentiable and non-zero 
on (0, 1), we see that u is twice differentiable in (0, 1) and hence a classical 
solution of equation (fTUI) . Considering the Frobenius expansions at the left- 
hand endpoint and the condition that u G L^((0, l),w{x)dx) we find that u G 
C°°([0, 1)) with u{0) = 0. Considering the Frobenius expansions at the right 
hand endpoint, we see that either u{x) ~ 1 as x ^ 1— or u{x) ~ (1 — x)~^^'^ 
as X — 1— . We are required to show that it is the former which holds. We 
have 

{Lf,u)^ = {f,L*u)^ (13) 



for all / G Dom(L). Since u is smooth in [0, 1), L*u 
the classical sense of differentiation. So 

— / {pf')'{x)u{x)dx = — f{x){pu')'{x)dx 
Jo Jo 



jJU 



-w 



-1 



{pu'Y in 



f {x)p{x)u' (x) 







= f'{x)p{x)u{x) — f {x)p{x)u' {x 

for all / G Dom (L). Since /(O) = n(0) = 0, this implies that 

' 



f'{x)p{x)g'{x)dx 

r-l 







lim p{x) f'{x)u{x) — f{x)u'{x) 



{pf'y{x)g{x)dx 



(14) 



for all / G Dom (A). By choosing f{x) = sin(7rx/2) we see that 

lim p(x)u'{x) = 0. 

X— >1 — 
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If u{x) ~ (1 — a;) ""^/^ as a; ^ 1— then u'{x) ~ (1 — x) ^ as x ^ 1— and 
hence p(x)u'(a;) ~ 1 as x — > 1—. This is a contradiction, so u{x) ~ 1 as 
required. 

We now assume < e < 1. It is immediate that (c) implies (d). The 
proof that (d) implies (a) is similar to the proof that (c) implies (a), but the 
possibility that u{x) ~ (1 — x)~^^^ as x — 1— is ruled out by the condition 
that M e L2((0, 1), w(x)dx). | 

Theorem 3.2 The Sturm- Liouville operator L is essentially self-adjoint. If 
< £ < 1 then Lc = L. 

Proof Suppose that fi is an eigenvalue of L*. Then, by Proposition 13. 11 /i is 
an eigenvalue of L and hence real, since L is symmetric. Hence the deficiency 
indices of L are both zero, so L is essentially self adjoint (see Theorem 1.2.7 
in [7]). If < e < 1 then L* is also essentially self-adjoint, by the same 
argument. Since L is a self-adjoint extension of Lc, the result follows. | 

Lemma 3.3 The operator L is injective. 

Proof Suppose for a contradiction that L is not injective. Then is an 
eigenvalue of L. By Proposition 13. H is also an eigenvalue of the classical 
Sturm-Liouville problem and hence of by the work in Section [21 
Davies showed in [2] that A > 1 for all real eigenvalues A of A^, so this is a 
contradiction. | 



4 Compactness of the resolvent 

In this section we give the integral kernel of the inverse of L explicitly, and 
use this to show that the resolvent is compact. This yields our result that 
the spectrum is discrete and the eigenvalues of L accumulate at -\-oo. 
We define 7 : [0, 1] ^ R U {00} by 

^(x) = r p{t)-^dt (15) 
Jo 

for all X e [0, 1] and G : [0, 1] x [0, 1] R U {00} by 

G{x,y) = \^f\ l^'ly (16) 
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Lemma 4.1 If G is as above we have: 

(i) G{x, y) eH for all x,y G [0, 1] except when x = y = 1; 
(a) G{x, y) = G{y, x) for all x,y G [0, 1]; 

(in) -§^G{x,y) = X[o,a)(l/)p(y)"^ for all x E [0, 1], y G (0, 1) \ {x}; 
(iv) £G{x,y) = X[o,y)ix)p{xy^ for all y e [0, 1], x G (0, 1) \ {y}. 

Proof (i) If X < 1 or ?/ < 1 then p^^ is bounded on [0, min{x, y}] and hence 
the integral is finite. 

(ii) Immediate from the symmetry of the definition. 

(iii) For y G (0,x), -^G{x,y) = i{y) = p{y)~^, whereas for y G (x, 1), 
lG{x,y) = ^^^{x) = 0.' 

(iv) Similar to the proof of (iii). | 

Theorem 4.2 The operator L has a compact inverse R given by 

(i?/)(x)= C G{x,y)f{y)w{y)dy (17) 
Jo 

for all f G L2((0, l),w{x)<lx) and all x G [0, 1). 

Proof We first prove that G G L^([0, 1] x [0,l],w(x)dx x w{y)dy), and 
hence that (fTTIl defines a Hilbert-Schmidt operator on L^((0, l),w{x)dx). If 
< y < X < 1 then 

\G{x,y)fw{y) = (^J^ Hv) 

< Co {{1 - yy'/' - if w{y) 

< Coy-' (l - (1 - y)'/') ((1 - yr'^' - l) (y + 1)'^/^ 
for some constant Cq. Hence 



\G{x,y)\'w{x)dxw{y) < coy-^ (l - (1 - y) V^) ((i _ y)-iA _ l) (y + ^(3;)dx 

< ciy-2 (l - (1 - y)'/') ((1 - y)-V- - l) (y + ly'/' f\l - x^'dx 



< C2y-' (l - (1 - y)'/') ((1 - y)-'/' -l){y + ly'^'il - y)'^''' 

< C2y-2(l_y)(^l_(l_y)l/e)'(y+l)-l/e 

for some constants Ci and C2 and all y G (0, 1]. As a function of y, this is 
continuous on (0, 1], and in a neighbourhood of we have 

\Gix,y)fw{x)dxw{y) < c,y-^/' ~ yiZ/f^^ ) ^ ^^y'^'^ 
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for some constants C3 and C4. Since e < 2 we conclude that 



1 n 2 

\G{x,y)\ w{x)dxw{y)dy < 00 

Jy 



and hence 



/ \G{x,y)\ w{x)dxw{y)dy = 2 / / \G{x,y)\ w{x)dxw{y)dy < 00 
JQ Jo Jo Jy 

as required, since G{y,x) = G{x,y) by Lemma WA\ 

We now prove that R is the inverse of L. Suppose that / G C^(0, 1). We 
have 

{Rf){x)= f' G{x,y)f{y)w{y)dy (18) 



and, since / is zero in sufficiently small neighbourhoods of and 1, it is easy 
to show that this is differentiable with 

= J^' ^G{x,y)f{y)w{y)dy 



X[o,y){x)p{x) ^f{y)w{y)dy 
Pixy^ I f{y)w{y)dy 



by Lemma 14.11 The last integral is smooth and vanishes in a neighbourhood 
of 1, and p{x)~^ is smooth on [0,1), so this implies that Rf G C°°([0, 1]). 
Also 

(i?/)(0)= (\{Q)f\y)w{y)dy = Q 



since 7(0) = 0. Therefore Rf G Dom (L) and 

— I p(x) — j 
dx \ dx Jo 



d / d ^ 
{LRf){x) = -w{xy^— [p{x)— G{x,y)f{y)w{y)dy 



-w{x) ^-^(pi^)J^ X[o,y)ix)p{x) ^f{y)w{y)dy 



= -w{x)-^— / f{y)w{y)dy 
dx Jx 

= /(^) 

for all / G C^{[0,1]). If / G V{{0,l),w{x)dx), let (/„) be a sequence in 
C,^([0,1]) such that ||/n-/IL as n ^ 0. Then \\Rf^-Rf\\^ ^ 
and \ \LRfn - /IL = ||/n - /IL as n ^ 00. Hence Rf G Dom (l) and 
LRf = f. 

Conversely, let / G Dom(L). Then RLf G Dom (l) and LRLf = Lf. 
Now RLf = f since L is injective by Lemma [3.31 | 
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Corollary 4.3 The Sturm- Liouville operator L is non-negative in the sense 
that Spec (jL^ C (0, oo) . 

Proof Since L has compact resolvent, it has empty essential spectrum, and 
since it is self-adjoint its spectrum is thus equal to the set of its eigenvalues. 
By Proposition 13.1^ it is sufficient to show that all eigenvalues of L are non- 
negative. If /i is an eigenvalue of L and / is a corresponding eigenvector with 
ll/L = lthen 

= (L/,/)^ = - [\pfJ{x)J{^dx = [\{x) \f'{x)\'dx > (19) 
JO Jo 

since p is non-negative on [0, 1]. Note that the inequality is strict, since /' = 
a.e. would imply that / is constant and hence 0, since /(O) = 0. | 

Corollary 4.4 There exists a complete orthonormal set of eigenvectors {/n}.^ 
of L with corresponding eigenvalues /in > which converge monotonically to 
+00 as n oo. 

Proof The corresponding result for L is standard, and the result for L follows 
by Proposition 13. 1[ | 

Corollary 4.5 The operator —iH defined in Section U\ has infinitely many 
eigenvalues which can be enumerated {^n}'^=-oo' increasing order, such 
that Xq = 0, \-n = — A„ and A„ ^ oo as n oo. 

Proof The eigenvalues of —iH are the same as the eigenvalues of A, since the 
two operators are unitarily equivalent. If {\n}'^=i are the eigenvalues of 
in increasing order then (jl]) tells us that the eigenvalues of A are {A„}^_g^, 
where Aq = and A_„ = — A„. It follows from Corollary 12.21 Proposition 13.11 
and Corollary 14.41 that A„ = efin/2 ^ oo as n ^ oo. | 

5 Eigenvalue asymptotics and numerics 
5.1 Quadratic form formulation 

When considering self-adjoint operators, it is standard practice to obtain 
eigenvalue asymptotics and upper and lower eigenvalue bounds by quadratic 
form techniques and variational methods. Our ffist task is to identify the 
precise domain of the quadratic form associated with L. We also show that 
C^(0, 1) is a form core for the associated form, which frequently allows us 
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to restrict our attention to this simpler class of functions in the subsequent 
analysis. 

We define a quadratic form Q on L^((0, l),w{x)dix) by 

W)= I'plffdx (20) 
Jo 

for / in the domain 

W'' = {/ e L^((0, 1), w(x)dx) : Qif) < 00} c {/ e C(0, 1) : /' e 
and define a norm 



\u\ 



{QH + \\u\\iy/' 



on 

Theorem 5.1 The domain W^'^ is complete with respect to ||-||^i,2. 

Proof Suppose that (/„) is a Cauchy sequence in W^'^. Let 6 e (0,1/2). 
Then 

"1-5 „ rl-S 



f 2 -1 f 2 

/ \fn{x)\ dx < sup w{y) / \fn{x)\ w{x)dx, 

JS y(Z^S^l-S) JS 

/ l/n(^)l < sup p(2/)"^ / \fn{x)\ P{x)dx 



y^iS,i-S) 

SO /n I (5,1-5) lies in the Sobolev space H^{6, 1 — 6). Similar inequalities show 
that (/n|(5,i-5)) is a Cauchy sequence in H^{5,1 — 5) and hence converges 
to some hs e H^{S, 1 — S). It is easy to show that if 61 < 82 then hs^ — 
/i5j(52,i-52)- Hence we may consistently define h,k : (0,1) — > C by h{x) — 
hi/m{x) and k{x) = h[^.^^^{x) for any m G N such that x G (1/m, 1 — 1/m). 
Clearly /;,, A; G -^^ioc(0! !)• Let G C^(0, 1) and choose m G N large enough 
that Supp (0) C (1/m, 1 — 1/m). Then 

rl rl—l/m 

/ /i(a;)0'(x)dx = / hi/m{x)(f)'{x)dx 

Jo Jl/m 

rl — l/m 

= - / h\i^{x)(t){x)dx 

Jl/m ' 

— — I k(x)(f)(x)dx. 
Jo 

Since this holds for aU G C^°°(0, 1), k = h'. 
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Since (/„) is Cauchy in W ' , Q{fn) and ||/n||^ are bounded. Let Ci 
sup„ Q{fn) and C2 = sup„ . For all m e N, 



l-l/m , 2 rl-l/m 2 2 

w(a;)da; < sup / \fn{x)\ w{x)dx < sup \\fn\\^ = C'a 



l/m 

and hence 
m e N, 

rl— l/m 
'1/: 



n Jl/ra 



< C2 by the Monotone Convergence Theorem. Also for all 



rl-l/m 2 rl-l/m ^ 

/ |/i p(a;)da; < sup / p(a;)da; < sup(5(/„) = Ci 

Jl/m n J\/m n 



and hence Q{h) < Ci by the Monotone Convergence Theorem. Therefore 
h e W^'"^. 

Let ?7 > be given. Since (/„) is a Cauchy sequence in W^''^, there exists 
e N such that 

/•I 2 

/ l/ni - /nal w(x)dx < 7^/6 
JO 

and 

for all ni,n2 > N. By the Dominated Convergence Theorem, 

/ \fN — h\ w{x)dx + / \fN — h\ w{x)dx ^ 

Jo il-l/m 

as m — > oo. Hence there exists M e N such that 

f^/m 2 2 

/ \fN — h\ w(x)dx + / l/jv — ^1 w(x)dx < 77/6 

Jo Jl-l/m 

for all m > M. li n > N then 

/„ — /iP t(;(x)da: < 2 / 



(0,l/m)U(l-l/m,l) 



'(0,l/m)U(l-l/m,l) 



|/n-/iv| W(x)dx 



+2 / |/jv-/irw(x)dx 

^(0,l/m)U(l-l/m,l) 
/•I 2 

< 2 / |/„ - /ivl w{x)dx 

Jo 

+2 / |/jv-/i|^'u;(a;)dx 

./(0,l/m)U(l-l/m,l) 

< 277/6 + 277/6 

= 277/3 
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for all m > M. Since w is bounded on (1/M, 1 — 1/M), 

l-l/M 



l-l/M 2 

\fn - h\ w{x)dx < sup w{y) 

1/M 2/e(l/M,l-l/M) 

= sup 'Ji'(y) 

2/e(l/M,l-l/M) 



as n — > oo. Hence there exists N' such that 

r-l-l/M 



1/M 
l-l/M 

1/M 



|/„-/i|'dx^O 



/n — ^1/M 



dx^O 



/ l/n - /i| w(x)dx < rj/3 

J 1/M 



for all n > N'. Therefore 



/ \fn - hf^ w{x)dx < T] 
Jo 



for all n > max(A^, A?"'), so — ^ as n — > oo. Similarly Q{fn — h) 



as n — > oo. Hence fn-^he W^'^ as n — > oo. ■ 

Theorem 5.2 If f E W^'"^ then f is continuous on [0, 1) and /(O) = 0. 

Proof For any a e (0,1), p-^^ e L2(0, a) and p^^V e L'^{0,a). Hence 
/' e iv^(0, a), so / is continuous on [0, a\. Since a is arbitrary in (0, 1), / is 
continuous on [0, 1). 

Suppose that /(O) ^ 0. Then \ ff > c on [0, 6] for some 5 > and some 
c> 0. Then 

/ w{x) \f{x)f dx > c w{x)dx — oo. 
Jo Jo 

However, 

['w{x) \f{x)\'dx^\\f 
Jo 

This is a contradiction, so /(O) = 0. | 
Theorem 5.3 The space C^(0, 1) is dense in W^''^. 
Proof For all 5 such that < 5 < 1/3, define fs by 



< oo. 



0<x<6 

f{2{x-6)) 6<x<2S 

f{x) 2S<x<l-S 

[ f{l-S) 1-S<x<l. 



13 



Then fs\{5,i) G 1) and fs{S) = 0. It follows from the boundedness of p 

and w on [6, 1] that fs e W^^'^. We have 

^\f{x)-fs{x)fw{x)dx = [' \f{x)fw{x)dx (21) 

Jo 

+ r\fix)-fi2ix-6))\'wix)dx (22) 
+ f \f{x)-f{l-5)\^w{x)dx (23) 

for all 5 G (0, 1/3). That the first integral on the right hand side converges 
to zero as (5 —i> is elementary. On some neighbourhood of 0, to is monotone 
decreasing. Thus, if we take 6 sufficiently small, then 'w{x/2 + 5) < w{x) for 
all X G (0, 26). Hence, for all small enough 6, 

2S „ r2S 1-25 

\f{x) - f{2{x-6))fw{x)dx < 2 \f{2{x-6))fw{x)dx + 2 \f{x)f w{x)dx 

J S J 6 

\f{x)fw{x/2 + S)dx + 2 / \f{x)fw{x)dx 

/•2<5 2 f"^^ 2 

< / |/(a;)| w{x)dx + 2 / \fx)\ w{x)dx 
Jo Js 



as 5 — s> 0. In order to with integral ( l23l) . we must first obtain an estimate on 
/ near 1. By Theorem 15.21 /(O) = 0. Hence 



|/(x)| < / !/'(.)! d. 

Jo 

< ^£p{sr'dsy' ^£p{s) \f'{s)fdsy 

for all X G (0, 1), so 

\f{x)f < Qif) rp{s)-'ds 

Jo 

< cQ{f){l ~ x)-'/^ 

for some constant c. Now, since w is monotone decreasing in a neighbourhood 
of 1 and there exists a constant c' such that w{s) < c'(l — s)-*^/"^ for all s > 1/2, 

\f{x) - f{l- 5)fw{x)dx < 2 \f{x)fw{x)dx 

l-S Jl-S 

+2 \f(l-5)fw(x)dx 

Jl-S 
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< 2 \f{x)fw{x)dx 

Jl-S 

+2 \f{l-5)fwil-5)dx 

Jl-5 

< 2 \f(x)fw{x)dx 

Jl-5 

+2cc'Q{f)6{l-6)-'/'{l-6)'/' 


as 6 0. Therefore || / — ^ as 5 ^ 0. Similarly, we may show that 

Qif — A) — as 5 — s> 0. Hence fs ^^—^ f as 6 ^ 0. 
Next we define fs^rj by 

r fs{x) 0<x<l-2r] 

A,(x)= - 2r/) (1 - l-2r^<x<l~rj 

[ l-r]<x<l 

for all S,r]< 1/3. Then E H\5, 1 - r/), /(x) = on (0,5) U (1 - 

?7, 1), and \i-m^^s+ f5,r){x) = lim^^i_^_ /^^^(a;) = 0. Hence fs,^ e if^(0,l). 
By Theorem 2 of Section 5.5 of |8j, there exists a sequence (/5,,,,n)neN in 

C^(5, 1 — 77) C C^(0, 1) such that fs,r],n — > /<5,r) as n 00. Since p and w 

are bounded on [5, 1 — r;], it follows that fs,r],n — ^ fs,ri as n ^ 00. 

It only remains to show that we may approximate fs by fs^n- For all 
(5 G (0,1/3), 

1/5(3;) - Ar, (a;) I ti'(x)dx = / \fs{x) - fs,rj{x)\ w{x)dx (24) 

+ [' \fsix)fwix)dx. (25) 
Clearly integral (125|) converges to zero as 77 — and 

\fs{x) - fs,r,{x)\ w{x)dx < 2 \fs{x)\ w{x)dx 

l-2r] Jl-2ri 

+2 /' ^ \fs,rj{x)f W{x)dx 

< 2 ['~^ \fs{x)fwix)dx 
+2 f' \fsil-2ri)fwix)dx 

< 2 f' \fsix)fwix)dx 

Jl-2ri 
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+2 





1-2J7 



w{x)dx 



as T] ^ 0. Hence 1 1/^ — fs,'q\ \ ^Oasr^^O. If < r] < 6 then 



/^(a;) - firjix) p{x)dx 



+ 

< 2 



+2 

< 2 



1— 2 
1 ^ 2 

1/5(0;) I p{x)dx 

1-V 

1--') 2 

p{x)dx 

l-2»7 



l-2r? 

1-'? 



'1-277 

+2r/-2|/,(l-2r^)r 



1-2J7 



p(x)da; 



< 2 / |/Kx)rp(x)da; 
+2^"MlMlioo sup 

l-2»7<z<l 

< 2 f'~'\n{x)fpix)dx 

Jl-2'q 

+2r^-^\\fs\\l^ (2r^)i+V.2i-Ve 
1-') 



l-2»7 





as ?7 ^ 0. Therefore /^^^ ^^^^ — > as 17 — > 0. | 

Corollary 5.4 The quadratic form associated with L is Q. 

Proof Let Q' be the form defined on Dom (L) by Q'{f, g) = {Lf, g)^. Then 
by Theorem 4.4.5 of [7], Q' is closable and its closure is associated with a self- 
adjoint extension of L, which must be L, since L is essentially-self adjoint. 
For / G Dom (L), 

QU) = (LfJ)^ 



{pfy{x)f{x)dx 
16 



= / p\f'\^dx 
Jo 

= Qif) 

so Q is an extension of Q. By the above theorem, 



C~(0, 1) C Dom [Q) C Dom (Q) = C^{0, 1) 
and hence Q is the closure of Q, as required. | 

5.2 Transformation to a Schrodinger Operator 

We next use a suitable change of variables to convert L into a Schrodinger 
operator on a certain space 1/^(0, /3). This allows us to use the extensive range 
of standard techniques available for controlling the eigenvalues of Schrodinger 
operators. 



Theorem 5.5 Define 

\ p{y) 



Jo 



dy 



for t e [0, 1], and let (3 = -0(1). Then p < oo and i/j : [0, 1] [0,/?] is 
invertible. Let (p — and define 

c{s) = iw{^{s))p{<P{s))r'/' 

v{s) = -<^)<^(^)^;ii;y-<^)<^(^)a-s^W 

Qc{9) = f^{\g'{s)\' + V{s)\g{s)f}ds 

for all g G C^(0,/3) C L^(0,/5). Then Qc is a closable quadratic form and 
its closure Q is associated with a self-adjoint operator H , which is unitarily 
equivalent to L. The potential V is smooth on (0, 1), V{s) ~ as s — > 0+ 

and V{s) ~ (i^ - |) (/3 - s)-^ as s ^ (3-. 



Proof We have 



/3 = Cy-"\i-yV\i + y)-"'dy 

Jo 

< V2 /'^' y-^/^dy + V2 [\l - yY^'^dy 
Jo Ji/i 

= 2^2 C'\-''^dy = A. 
Jo 



17 



Since w{t)^^'^p{t)~^^'^ is smooth and positive for all t G (0, 1), it is immediate 
from its definition that ip is smooth on (0, 1) and continuous and strictly 
monotone increasing on [0,1]. Thus ip is injective. Since ip is continuous, 
ip{0) = and = (3, ip must be surjective. Hence -0 is invertible. It is 
easy to show that (j) and (j)' are smooth and non-zero. 
We define 

{Uf){s) = c{s)-'f{<P{s)) 

for all / G L^((0, l)/u}{x)dx) and all s G [0, 13]. It follows from the definitions 
of tP and that i^^t) = w{t)^/^p{t)-'^/'^ and (f)'{s) = p((/.(s))i/2^((/.(s))-V2. 
Hence, making the change of variables x — (f){s), 



2 m2 



w 







V(o) 





Jo 

= l^\c{s)-'f{cl>{s))\'ds 



for all / G L^((0, 1), u'(x)da;). It follows that t/ is a unitary operator from 
L2((0,l),w(x)dx) to L2(0,/3). 

Since p and w are smooth and non-zero on (0, 1) and is smooth and 
non-zero on (0,/9), c is smooth and non-zero on (0,/3). Hence Uf G C^(0,/3) 
if and only if / G C^(0, 1). Making the change of variables x — (j){s) as 
above, 

Q{f) = ['\f{x)\'p{x)dx 
Jo 

= [^\c\s){Uf){s) + c{s){Ufy{s)U'{cl>{s))MHs))cj>\s)ds 
Jo 

^\c'{s){Uf)is) + cis)iUfyis)U'is)-Wis))ds 

( [c{sf \{Uf)'{st + c'(.)c(.)^ {\{Uf){s)f) + c'{sf ^^d. 
£ {c{sf\{Uf)'{s)\' + (c'(.)^ - {cc'y{s)) \{Uf){s)\'] ^ 



p(0(s)) 



<P'{s) 
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'{mfns)\' + vis)\{uf)is)\'}ds 

= QciUf) 

for all / eC^{0,l). 

Since Qc is the form arising from the symmetric operator 

H^ = -A + V 

with domain Dom (ifc) = C^(0,/5) and 

{HJJ)=Q,{f) = Q{U-'f)>0 

for all / G C^(0,/3), Qc is closable and its closure, Q, is associated with a 
non-negative self-adjoint extension i7 of He by Theorem 4.4.5 of [7]. 
We have proven that / G C;?°(0, 1) if and only if Uf G C^{0, p), and 

= QciUf) (26) 



for all / G C^(0, 1). We now prove that f// G Dom (^Q j if and only if 
/ G Dom (Q) and that 

= Q{Uf) (27) 

for all / G Dom((5). By Theorem 15.31 C^{0, 1) is dense in Dom (Q) with 
respect to |Hlvi/i,2. Let / G Dom{Q). Then there is a sequence (/„) in 

C,^(0, 1) such that fn-^f as n ^ oo. Hence (/„) is a Cauchy sequence in 
C^(0, 1) with respect to ||-||y(^i,2. By f l26l) and the unitarity of U, {U fn) is a 
Cauchy sequence in C^(0, /5) with respect to |||-||| = (|Hli2 +<5("))"'^^^- Since 
Dom (^Q) is the completion of C^(0,/5) with respect to |||-|||, there exists 
g G Dom i^Q^ such that — (^lll — * as n ^ oo. Also H/n — / ||^ — ^ 

as n — oo implies that ||f//„ — f//||^2 ^ as n oo. It follows that 
Uf = (7 G Dom (q)- The converse is similar. For all / G Dom{Q), 

Q(uf) = \\\uf\\\'-\\uf\\l, 

= hm WWUW?- hm ||f//„||i2 

n— >oo n— »co 



|2 llxl|2 



= Q(/) 

where, as before, (/„) is a sequence in C,°°(0,1) such that /„ f. 
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It now follows from the polarisation identity for sesquilinear forms that 

Q'{Uf,Ug) = Q\f,g) 

for all f,g & Dom{Q). It follows immediately that Uf G Dom (if) if and 
only if / G Dom (^L^, and that H = ULU^^. Since c, c', c", p, (j) and 1/0' are 
smooth on (0, 1), it follows that V is smooth on (0, 1). For the asymptotics 
of V , see Appendix [X] | 

5.3 Eigenvalue asymptotics 

Throughout this section, {iJin)^=i shall be the eigenvalues of L, or equivalently 
of L by Proposition listed in increasing order and repeated according to 
multiplicity as in Corollary 14.41 By Theorem 15.51 L is unitarily equivalent to 
a Schrodinger operator H with Dirichlet boundary conditions. We shall use 
the Rayleigh-Ritz variational formula to obtain bounds on /x„ in terms of the 
Dirichlet eigenvalues of —A on various intervals. Recall that the Dirichlet 
eigenvalues of —A on the interval [a, 6] are {n^7r^(6 — a)^^ : n G N} with 
corresponding eigenf unctions {sin(n7r(x — a)/{b — a)) : n G N}. We quote 
the variational formula from [7J: 

If i^" is a non-negative self-adjoint operator on a Hilbert space Ti and M 
is a finite-dimensional subspace of Dom (K) then we define 

A(M) = snp{{Kf,f) : / G M and 11/11 = 1} 

and 

A„ = inf{A(M) : M C Dom (K) and dim(M) = n}. 

If A„ —* +00 as n — s> oo then K has compact resolvent and the numbers A„ 
coincide with the eigenvalues of K written in increasing order and repeated 
according to multiplicity. Since K is non-negative and self-adjoint, it is 
associated with a closed quadratic form Q. If P is a core for Q, that is, a 
subspace of the domain of Q such that the closure of Q restricted to P is Q, 
then we have 

A„ = inf{A(M) : M C D and dim(M) = n} 

for all n G N. 

Theorem 5.6 Let f3 = /g^y"-'^/^(l —y)^'^^'^{l + yy^^'^dy, as in Theorem \5.5[ 
Then ^ 

lim n"V„ = — . 
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Indeed, if a = mm{V{s) : s G (0,/3)}, where V is as in Theorem \5.5l then 

[In > + a 

for all n G N, and 

as n ^ oo. 
Proof Define 

Qa{g) = {Hg, g)-a {g, g) = Qdg) -a{g,g)>0 

for all g G C^(0,/5). Then Qa is a closable form and its closure is the 
form associated with a non-negative self-adjoint extension Ha of He — al, 
where H^ = -A + V with Dom (He) = C^{0, Since C^{0, p) is a core for 
both Q and Qa, the variational formula implies that the nth eigenvalue of 

Ha is Hn — OL. 

We now define 

Q{g)= f\g\s)f^s< f \g\s)f ^ iy^s) - a)\g{s)\Us = Qa{g) 
Jo Jo 

for all g G C^(0,/9). The closure Q of Q is the form associated with the 
operator —A on L^(0, (3) with Dirichlet boundary conditions. Since C^(0, (3) 
is a core for both Q and Qa, the variational formula implies that 

^ < /in - a 

for all n G N. 

For all sufficiently small S > define 

Ks{f) = {-A + csI){f) 

for all / G C,°°(5, /5 - 5), where 

Cs = sup{l^(s) : s E {6, P — 6). 

Then Ks is a non-negative symmetric operator and hence has a non-negative 
self-ajoint extension Kg, called the Friedrichs extension of K^, associated 
with the closure of the form Qs defined by 

Qs{f) = {KsfJ) 
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for all / e C^{S, (3 — 6). For all finite dimensional subspaces M of Dom {K^), 
define 

\s{M)^SM-p{{Ksf,f) : / e M and ll/ll = 1} 

and 

\n = inf{A5(M) : M C C~(5, (3-5) and dim(M) = n}. 
Then, for MCC,°°(5,/3-5), 

A5(M) = sup{((-A + c,/)/,/):/GMand 11/11 = 1} 
= sup{((-A/,/):/eMand 11/11 = 1} + C5 



and hence 



by the variational formula for the Dirichlct eigenvalues of —A on {6,(3 — 5). 
Since A^^^ — > +oo as n — > cxd and /3 — 5) is a core for the quadratic form 
associated with Ks, Ks has compact resolvent and its eigenvalues coincide 
with Xs,n- For all M C Dom (H), define 

/x(M) = sup{{Hf,f) : / e M and ll/ll = 1}. 

li f eC^(5, (3 -5), then 

{Hf, /) = ((-A + V)f, f) < ((-A + csl)f, f) = {Ksf, f) 

so /i(M) < Xs{M) for all M C C^{S, (3-5). Since C^°°(0, /3) is a core for Q, 

tin = inf{/i(M) : M C C,°^(0, /?) and dim(iV/) = n} 

< inf{^(M) : M C C^{S, (3 - S) and dim(M) = n} 

< inf{A5(M) : M C C^{6, l3 - 6) and dim(M) = n} 



+ Cs 



for all n e N. From the asymptotics of V, csS'^ max{|, ^ — =: c as 
5 — > 0. Hence, given u > 0, cs < {c + v)5~'^ for all sufficiently small 5. For 
such 5, 

^^^^ 22/^ 1 l\ C+t/ ^/c-\ 
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for all n G N. Hence, for all sufficiently large G N, 



/52 



2 2 4n-2/3(/5-n"2/3) 



/?3 



as n ^ oo. 



5.4 Eigenvalue numerics 

Let {a*^™)}5^^]^ be a monotone decreasing sequence in (0, 1) converging to 
and {h^'^^}'^^i a monotone increasing sequence in (0, 1) converging to 1 such 
that a^^^ < b^^\ For each m G N, we consider the operator L*^™"^ defined by 
ffT^ on the domain 

= {f e Dom (l) : G C°°(/("^)) and /(a^™)) = /(fe^'")) = 0}. 

Then L^"^^ is the operator corresponding to a regular Sturm-Liouville prob- 
lem. By Theorem 14 of [9], L^™'^ has infinitely many eigenvalues /x^'"-' tending 

to +00. 

Lemma 5.7 For each n G N, {/ii'^^}m=i ^■^ decreasing sequence such that 
f^^rT^ > fJ'n for a// m G N. 

Proof It is easy to see that 

c r'^^) c . . . C Dom (Z) 
and hence the result follows immediately from the variational formula. | 
Lemma 5.8 For each n G N, fi^^ — > yUn as m ^ oo. 
Proof Since L has compact resolvent and C^(0, 1) is a form core for L, 
Hn = inf {/i(M) : M C C,°°(0, 1) and dim(M) = n} 

where /i(M) is as in the variational formula. Hence, given rj > 0, there exists 
M C C^{0, 1) such that dim(M) = n and /i„ < /i(M) < /i„ + r/. Since M is 
finite-dimensional, we can choose mo G N such that Supp (/) G (a^™"-*, ^('"o)) 
for all / G M. Hence M C T'('^o) and so /i^T*") < /i(M) < fi^ + V- Combining 
this with Lemma 15.71 we have yU„ < /i^*"-* < fin + for all m > mo . 
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Theorem 5.9 If Xn is as in Corollary \4.4\ and A^™'^ = e/i^^V^; then, for 
each n e'N, decreases monotonically to Xn as m oo. 

Proof Immediate from Lemmas 15.71 and 15.81 | 

We have shown that we can obtain upper bounds on the eigenvalues of 
L and —iH by computing the eigenvalues of regular Sturm-Liouville prob- 
lems, and moreover, that these bounds converge to the eigenvalues of these 
operators as we allow the endpoints of the regular problems converge to the 
endpoints of the Sturm-Liouville problem associated with L. We give ta- 
bles of numerical calculations, obtained using the software package SLEDGE 
[To] (obtainable from|http:/ /www.netlib.org/misc/sledge ), of A^™^ for a^"^) = 



10"™, 6^"^^ = 1 — 10"™ and various values of e. For comparison, we also 
include numerical calculations from [T], [2] and [3]. Our computations are 
performed with an absolute error tolerance of 10~^, so we can be reason- 
ably confident that any calculated values which exceed our upper bounds by 
more than this amount are not accurate to the stated precision. All values 
given in [2] are consistent with these bounds when rounded to two decimal 
places, but there are some discrepancies at higher levels of precision. All 
values given in [T] and [5] are consistent with these bounds when rounded to 
one decimal place, but both sets of calculations have discrepancies at higher 
levels of precision. 

A Asymptotics of V 

We want to analyse the asymptotic behaviour of 

V(s) = -c(sy(sf-^ _ ,,,),,(,)£(£M^(f|_fW^ (28) 

as s — > 0-|- and s — > /3— . We have 

c = fc|(o,i)O0 (29) 
c' = (A;'|(o,i)O0)0' (30) 
c" = (A;"|(o,i)O0)0'2 + (A;'|(o,i)O0)0" (31) 

where k is the analytic function defined by 

k{z) = Z^l\\ - z)-^'^'-^/\l + z)V2.-V4 

for all z e C \ ((-oo, 0] U [1, +oo)). 
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n 


;\(3) 

n 


A„ 


;\(5) 

n 


;\(6) 

n 


n 


A„ [5] 


A„ [3] 


1 


1.45457 


1.44906 


1.44851 


1.44845 


1.44844 


1.4485 


1.449323 


2 


4.34574 


4.31891 


4.31614 


4.31587 


4.31584 


4.3159 


4.319645 


3 


8.70318 


8.63035 


8.62264 


8.62186 


8.62178 


8.6219 


8.631474 


4 


14.53324 


14.38251 


14.36590 


14.36421 


14.36405 


14.3638 


14.382886 


5 


21.84048 


21.57464 


21.54473 


21.54167 


21.54137 


21.5414 





Table 1: Numerical calculations for £ = 1, compared with calculations from 
[2] and [3j 



n 


A(3) 


A(4) 


A(5) 


A(6) 


A(7) 


An [3] 


1 


1.17382 


1.16782 


1.16720 


1.16714 


1.16714 


1.167342 


2 


2.99250 


2.97016 


2.96847 


2.96823 


2.96821 


2.968852 


3 


5.54084 


5.48803 


5.48231 


5.48174 


5.48168 


5.483680 


4 


8.82509 


8.72519 


8.71398 


8.71284 


8.71272 


8.715534 


5 


12.85050 


12.68265 


12.66336 


12.66138 


12.66119 




6 


17.61828 


17.36052 


17.32987 


17.32674 


17.32643 




7 


23.13086 


22.75976 


22.71552 


22.71081 


22.71033 




8 


29.39064 


28.88240 


28.81847 


28.81174 


28.81106 




9 


36.39780 


35.72664 


35.63949 


35.63022 


35.62928 




10 


44.15374 


43.29376 


43.17838 


43.16790 


43.16666 





Table 2: Numerical calculations for e = 0.5, compared with calculations from 



n 


A(3) 


A(4) 


A^(f) 


A(6) 


A(7) 


A„ \V\ (numerical) 


A„ 


1 


1.02908 


1.01149 


1.00961 


1.00942 


1.00940 


1.0097 


1.00968 


2 


2.11378 


2.07759 


2.07349 


2.07306 


2.07305 


2.0733 


2.07334 


3 


3.29583 


3.23676 


3.22974 


3.22902 


3.22894 


3.2297 


3.22978 


4 


4.59835 


4.51260 


4.50208 


4.50099 


4.50088 


4.5012 


4.50134 


5 


6.03392 


5.91589 


5.90082 


5.89984 


5.89968 


5.8992 


5.89993 


6 


7.60918 


7.45354 


7.43391 


7.43175 


7.43154 


7.4298 


7.43194 


7 


9.32789 


9.13017 


9.10350 


9.10063 


9.10034 


9.0951 


9.10097 


8 


11.19231 


10.94654 


10.91287 


10.90919 


10.90881 


10.8945 


10.9092 


9 


13.20382 


12.90464 


12.86256 


12.85789 


12.85742 


12.8252 


12.8578 


10 


15.36360 


15.00536 


14.95367 


14.94786 


14.94727 


14.8820 


14.9478 



Table 3: Numerical calculations for e = 0.1, compared with calculations from 
[I] and [2] 
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Lemma A.l Asymptotically, 



1/4 



as z ^ and 



k{z) ~ z 
k'iz) ~ -z-^l^ 



-1/4 



k"{z) ~ + ^) + ^) 2^/^^-^/'(l - ^)-^/^^-^/^ 

as 2 — s> 1. 

Proof We may write 

k{z) = z'^^koiz) 

where 

ko{z) = (1 - ^)-i/2--i/4(l + ^)i/2--i/4 (32) 

is analytic in a neighbourhood of 0. The function k^ has a power series 
expansion 

oo 

^o(^) = 

n=0 

vahd in some disc -D(0; i?o)- Putting 2; = in flH^ . we see that Oq = 1. We 
thus have 

00 

1/4 ^ \- a„z"+i/^ 



fc(^) = .^/^ + E 



k\z) = lz-'/' + ^{n + l) 



n=l 

00 

n-3/4 



Qjrt.Z 



k"{z) = -^^,-V^ + ±{n + ^^[n-\]a.z--^ 

in the cut disc D(0;i?o) \ (— -Ro^O], and the stated asymptotics as 2; — > 
follow. 

Similarly, by writing 

k{z) = {i-zy^/^'-^'^ki{z) 



26 



where 

is analytic in a neighbourhood of 1, we obtain the stated asymptotics as 



Lemma A. 2 Asymptotically, 



i)'{s) ~ 2-\s 



as s 0+ and 



s 



\2 



<P"{s) ~ -1 

as s — 13—. 

Proof The inverse ip oi (f) is given by 
for all t e [0,1]. Since 







w{yY/My)-'^' = y-'^\i - y)-'^\i + yV 

for all y G (0,1), 
as t 0+ and 



as t ^ 1—. Putting t = 4>{s), we obtain the stated asymptotics for cj). 

The stated asymptotics for follow immediately from the asymptotics 
of (f) and the fact that 

(33) 



^'(0(s)) 

= pi<l>is)y/MHs))-'^' (34) 

= 0(.)V2(l_^(,))l/2(i + ^(,))l/2 (35) 

for all s e (0,/3). Differentiating (1351) and using the asymptotics we have 
calculated for and 0', we obtain the stated asymptotics for 0". | 
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Corollary A. 3 Asymptotically, 

c{s) ~ 2-'/h'/^ 

C'{S) ~ 2-3/2^-1/2 

c"(s) 2-^/2^-3/2 



as s — > 0+ and 



c{s) ~ 2i/^(/3-s)-i/^-i/2 



c"is) 



1 1^ 

1 1\ /I 3^ 



c'{s) ~ -+ 2i/^(/5-s: 



as s — >• /3— . 

Proof This follows immediately from equations (l29i) . ( l30l) . ( l3Ti) and Lemmas 
EUandlD I 

Corollary A. 4 Asymptotically, 

p{(f){s)) ~ 1 

p'(0(s)) ~ 

as s — * 0+ anc? 

y(0(s)) ~ -{i+^y'-'/^p-sf/^ 

as s ^ (3—. 

Proof The asymptotics of p and p' are: 

p{x) ~ 1 
p (a;) ~ — 

as a; — >• and 

p{x) ~ 2^-1/^(1 -x)i+i/^ 
p'(a;) ~ - + 2^-1/^(1 -x)i+i/^ 



as X ^ 1. The stated asymptotics follow from this and Lemma [A. 2[ | 
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Theorem A. 5 Asymptotically, 




as s ^ 0+ and 
as s ^ 13—. 

Proof This foUows from fl2Hl) . Lemma rA.2l and CoroUaries IA.3I and IA.4I | 
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